Phase Control of Nonadiabaticity-induced Quantum Chaos in An Optical Lattice 
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The qualitative nature (i.e. integrable vs. chaotic) of the translational dynamics of a three-level 
atom in an optical lattice is shown to be controllable by varying the relative laser phase of two 
standing wave lasers. Control is explained in terms of the nonadiabatic transition between optical 
potentials and the corresponding regular to chaotic transition in mixed classical-quantum dynamics. 
The results are of interest to both areas of coherent control and quantum chaos. 
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Recent years have witnessed an increasing interest in 
the coherent control of atomic and molecular processes 
. One central aspect of coherent control is phase con- 
trol, in which different optical phases are introduced into 
coherent laser-matter interactions in order to manipulate 
quantum interference effects and thus to achieve target 
objectives. It has been shown that phase control ap- 
proaches are widely applicable ]jj , even to some systems 
displaying quantum chaotic dynamics 0. 

As we show in this letter, optical lattices, of great re- 
cent interest 0] , provide an important system in which to 
explore aspects of quantum chaos and coherent control. 
In particular, we describe the all-optical phase control 
of the translational motion of atoms in a system that 
allows chaotic vs. integrable motion depending on the 
phase shift between two standing-wave laser fields. The 
control mechanism is shown to originate in the nona- 
diabatic coupling between different optical potentials as 
well as in the regular to chaotic transition in a mixed 
classical-quantum description of the model system. The 
results are of broad interest to both coherent control and 
quantum chaos. 

Consider a A-type 3-level atom moving along two co- 
propagating standing- wave laser beams, with two lower 
degenerate levels |1) and |3), and one upper level |2). 
Two laser fields, with different polarizations er+ and cr_, 
couple |1) with \2) and \2) with |3), respectively. A closed 
3-level A configuration of this kind may be realized, for 
example, in 4 He using the 2 3 Si — > 2 3 Pi transition ||. 
The laser fields are of the same frequency, with large de- 
tuning A from |2). We use x, p, M, Qi (f2 2 ), k\ (= k 2 ) to 
represent the position, momentum, atomic mass, the two 
Rabi frequencies (assumed real) and the two wavevec- 
tors, respectively. The relative phase of the two standing- 
waves is denoted by <f>. For generality we employ a set of 
natural units by scaling all the parameters, i.e., x° = A 
for x, p° = h/\ for p, t° = MX 2 /h for the time variable t, 
il° = 2/t° for A, Sli and f^- In terms of these units, the 
dynamical equations do not explicitly contain the atomic 
mass, the effective wavevector is given by k = 2tt, and 
[x,p] = i in the full quantum dynamics. In the rotating 
wave approximation and in the interaction picture, the 



Hamiltonian describing the translational motion along 
the laser beams is H = p 2 /2 + (Vij), with the potential 
matrix (Vij) (i,j = 1,2,3) given by 
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FIG. 1. Time dependence of (p) obtained from quantum 
wavepacket dynamics calculations (solid lines) and the ensem- 
ble statistics in a mixed classical-quantum description (dashed 
lines). The initial internal state is |1), and the (x) and (p) of 
the initial Gaussian ensemble are 0.0 and 25.0, respectively. 
The initial variances in position and momentum are chosen 
to be 1.0/10\/2 and 10.0/\/2, respectively. The relative phase 
cj> equals 0.0, 0.25-7T, and 0.57T in (a), (b) and (c), respectively. 

To demonstrate phase control of the dynamics, con- 
sider numerical results for the specific case where (1) 



1 



fii = 6.0 x 10 3 , f2 2 = 7.0 x 10 3 , A = 1.5 x 10 4 ; (2) 
the atom is initially in internal state |1); and (3) the av- 
erage momentum, average position, momentum variance, 
and position variance of the initial Gaussian wavepacket 
are given by (p) = 25.0, (x) = 0.0, dp = 10/V2 and 
Sx = I.O/IOa/2, respectively. In 4 He, this corresponds to 
t° ~ 77 /usee, the detuning ~ 27T-62 MHz (about 38 times 
the linewidth of |2)), and the initial kinetic temperature 
~ 29 [xK. Results for this case represent typical observa- 
tions for a wide range of system parameters and initial 
conditions that we have examined || . 

The solid curves in Fig. [ij which contain the essen- 
tial result of this letter, show the time dependence of 
the quantum momentum expectation value (p) for time 
evolving wavepackets for various <j>. The <j> = case [Fig. 
[j]a] displays a perfectly regular recurrence pattern. By 
contrast, the <j) — 0.25-7T case [Fig. [j]b] is characteristic 
of a relaxation process, with significantly irregular os- 
cillations of small amplitude. In addition, the associated 
power spectrum (not shown) is quite noisy, characteristic 
of chaos . Further tuning <j> leads to totally different dy- 
namics: in the <fi = 0.5-7T case [Fig. |]c], (p) lies very close 
to its initial value and undergoes regular oscillations, but 
with a characteristic frequency that is much higher than 
that in the <j> = case (The dashed curves in Fig. |l| 
are discussed later below). Clearly, the atom's transla- 
tional motion undergoes significant qualitative changes 
with controlled changes in <f>. 
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FIG. 2. The sensitivity of the quantum dynamics to slight 
changes of the relative phase parameter cf>. \ is the absolute 
value of the overlap between the two time evolving wavefunc- 
tions emanating from the same initial state as in Fig. [jj with 
the relative phase of the two standing-wave laser fields given 
by <j) and <f) + n/400. 



To demonstrate that this is indeed an integrable to 
chaotic transition, we consider the sensitivity of the dy- 
namics to slight changes of <p. Figure |^ shows the time 
dependence of the absolute value of the overlap \ of two 
time evolving wavefunctions emanating from the same 
initial Gaussian wavepacket, with the relative phase of 
the two laser fields given by <j> an d 4> + 7r/400, respec- 



tively. For (f> = or <f> = 0.57T, \ remains near unity 
throughout, indicating that the dynamics is insensitive 
to tiny changes of <j). By contrast, for <fi = 0.25-7T, \ is 
already less than 0.70 at t = 0.7. This interesting hy- 
persensitivity j§] to perturbations in the relative phase 
parameter 4> further suggests that the 4> — 0.25-7T case is 
indicative of quantum chaos. 
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FIG. 3. Time dependence of momentum for classical tra- 
jectories obtained by solving Eqs. (0), (^), and (Q), with the 
initial position x = and the initial momentum p = 25.0. 
The initial internal state is given by 1 1) . <f> equals 0.0 , 0.257T, 
and 0.5-7T in (a), (b) and (c), respectively. 

To further substantiate that this is a regular to chaotic 
transition we consider a mixed classical-quantum de- 
scription of the dynamics, i.e., where the center-of-mass 
motion evolves classically on an average potential and 
where the internal motion is treated quantum mechani- 
cally Jq|,P~0[| . From this perspective, the state of motion is 
described by a phase space point (x,p) and an internal 
wavefunction \ip) = (CijC^Ca), where Ci, C2, and C3 
are the projections of the wavefunction onto the three 
internal levels. The dynamical equations are given by 
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-dt =P ' 
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= > VkiCi, k = 1, 2, 3. 
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One can readily solve these equations numerically to 
obtain for the classical translational motion. The 
results, for the same system parameters and the same 
initial internal state as in the quantum calculations, but 
for a single trajectory initially at (x,p) — (0.0,25.0), are 
shown in Fig. ^. In particular, for (ft = [Fig. |]a] the 
oscillation of momentum is perfectly regular, in agree- 
ment with the regular recurrence pattern in Fig. [TJa. By 
contrast, for <ft = 0.257T [Fig. ||b], the trajectory is highly 
irregular, with random alternations between fast small- 
amplitude and slow large- amplitude oscillations. In this 
case two initially nearby trajectories shows exponential 
divergence in phase space with an associated Lyapunov 
exponent of ~ 50/t°. This supports the view that the 
irregular dynamics in Fig. and the sensitivity of the 
dynamics to slight changes of (ft shown in Fig. [^, are in- 
deed due to optical-phase-assisted quantum chaos. Fur- 
ther, in the (ft — 0.57T case [Fig. ^p], the regular classical 
motion is restored, with a characteristic frequency iden- 
tical to that in Fig. |l|c. Hence optical-phase control is 
evident in this classical-quantum treatment as well. 

An ensemble statistics in the classical treatment of 
translational motion provides further support. The 
dashed lines in Fig. |l| display the time dependence of 
the average momentum (p) for an ensemble of trajecto- 
ries initially centered at x — and p — 25.0, with the 
same initial variances as in the quantum calculations. 
Each individual trajectory in the ensemble is obtained 
by solving Eqs. @ - (^). The quantum-classical corre- 
spondence for the regular dynamics at times t < 0.15 in 
Figs. |l|a and [ljcj is impressive. On the other hand, as seen 
in Fig. |j]b, (p) for the classical ensemble quickly relaxes 
to zero, whereas (p) for the quantum ensemble remains 
far away from zero. This quantum-classical difference 
constitutes an excellent example of quantum suppression 
of classical chaos in an unbounded Hamiltonian system. 

Insight into the origin of phase control can be obtained 
by considering the dynamics in an adiabatic representa- 
tion. To do so we introduce an orthogonal transformation 
(Oij) (i,j — 1,2,3) to diagonalize the potential matrix 
(Vij). (Oij) is given by 
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where £(x, 



SI 2 sin 2 (kx) + Q 2 sm2 (kx + (ft) and 



r)(x, (ft)= ■J £ 2 (x, (ft) + A 2 . Corresponding to the three 
eigenvectors (Oij, 02j, O^j) (j = 1, 2, 3), are three eigen- 
potentials Vi(x,(ft) given by V\(x,(ft) = A + rj(x,(ft), the 
constant potential V% = 2 A, and Vz(x, (ft) = A — r](x, (ft). 
For the special case of (ft — 0, V\ and V% are degener- 
= nir, and the eigenvector (012,022,032) is 
x-independent; for the general cases of (ft ^ 0, V\(x, (ft) > 
V-2 > Vz(x, (ft), i.e., the three potential curves do not cross 
one another. Of particular interest is the constant poten- 
tial V2, associated with the eigenvector (012,022,032). 
Since 



Sli sin(fca;)Oi2 + ^2 sin(fcx + (ft)0. 
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V2 results from the complete quantum destructive inter- 
ference between the two standing-wave laser fields. As 
such, V% is an extension of the "dark optical lattice" in 
the presence of two counter-propagating plane- wave laser 
beams [Ulj,|lJ. 

Consider now transforming Eqs. (||) and (Q) to the 
eigen-potential (adiabatic) representation. Specifically, 
consider the dynamics in terms of Ci, i = 1,2,3, where 
&i = Yl=\OkiC k - Using J2l=i OkiOkj = Sij and 
OkidOkj/dx = ~X)fe=i OkjdOki/dx, Eqs. (||) and 

) can be transformed to 
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where the diagonal terms are the three adiabatic poten- 
tials given above, and the potential coupling terms are 
given by h 2 (x,(ft) = pkfl^ sm((ft) / [2^ 2 (rj 2 - A77)] 1 / 2 , 
t 13 (x, (ft) = pkA[fll sm(2kx)+nl sin(2fca;+20)]/4n 2 ^, and 
ha(x,<f>) = pktt^ sm((ft)/[2e(ri 2 + &v)} 1/2 A Note 
that tij (i =/= j = 1,2,3) is proportional to the momen- 
tum p. Thus, the coupling between optical potentials is 
due to the nonadiabatic effects associated with transla- 
tional motion. Further, for the computational example 
discussed above (and many other cases in which A > 0) 
one has mm[Vi(x,(ft) — V2] « min[V2 — V 3 (x,(ft)} and 
ti 2 (x, (ft) » t 2 3(x,(ft), suggesting that Va(x, (ft) is effec- 
tively decoupled from Vi (x, (ft) and V2. Note also that at 
the initial location x = 0, the initial internal state |1) is 
a superposition state of the two eigenvectors associated 
with V\ and V2 for (ft — 0, and reduces to the eigenvector 
associated with V2 for eft ^ 0. 

The key role of the relative laser phase (ft becomes clear 
as one compares the magnitude of the nonadiabatic cou- 
pling term ti 2 (x,4>) with that of (Vi — V2). For case 
(a), (ft = and ti2(x, (ft = 0) = 0, so the dynamics 
is adiabatic. Thus, in the case of Fig. [l]a, the quan- 
tum ensemble divides into two sub-ensembles: one is 
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trapped in one well of V% around x = and under- 
goes periodic oscillations, and the other experiences the 
trivial motion on the constant potential V%. For case 
(b), <j> = 0.257T. Here the smallest gap between V\ 
and V% is given by g(<fr) = \J A 2 ((f>) + A 2 — A, where 
A 2 ((f)) = [ft 2 +n| - y/tlf + Of + 2Sl 2 ft 2 cos(2</>)]/2. The 
corresponding ratio of the potential coupling term t\2 to 
g(<t>) is given by 



T(0) = 



pkQi&,2 sin(</>) 



V2[v4 2 (</>) + A 2 ] 1/4 ^(<£)ff 3/2 (0) • 



(9) 



Since T(</> = 0.257r) 1.0, the magnitude of the nonadia- 
batic coupling is comparable to that of (V1—V2), resulting 
in strong nonadiabatic effects. Thus, the chaotic motion 
seen in Fig. |3|b is induced by the significant nonadia- 
batic coupling between the two simple one-dimensional 
potentials V\ and V2. Finally, for case (c), (f> = O.bir. 
Here T(</> = 0.57r) = mm[T(<fi)} ss 0.16, i.e., the nonadia- 
batic coupling is appreciably weaker than in the case of 
(f) = 0.257T. As such, the translational motion, initially 
launched on the adiabatic potential V% , would essentially 
remain on V2 , with small perturbations from the insignif- 
icant Rabi population oscillation between Vi and V2. To 
further confirm this picture, one finds that the charac- 
teristic frequency of the regular dynamics in Fig. [lp and 
Fig. is ~ 1425, a value consistent with the Rabi fre- 
quencies given by \J (V\ — V2) 2 + 4i 2 2 [see Eq. (||)]. 

A number of additional remarks are in order. First, 
the two-standing-wave configuration is essential in this 
system. That is, if either or both of the two standing- 
wave fields are replaced by a traveling- wave, the nonadi- 
abatic coupling or the spacing between V\ and V2 is no 
longer a sensitive function of <fr, and there is no significant 
phase control. On the other hand, a three-level atom in 
two standing-wave laser fields of different but commen- 
surate frequencies also shows dynamics which is control- 
lable by changing cf> [Q. Second, in contrast to some 
recent studies on optical- magneto lattice , the cou- 

pling between different optical potentials discussed above 
is not due to additional magnetic fields, but directly 
due to nonadiabaticity. Further, unlike the work in Ref. 
fbSf , here we have observed clear signatures of quantum 
chaos in the quantum dynamics. Thus, this model is 
the first all-optics realization of nonadiabaticity-induced 
quantum chaos, a phenomenon first discovered in molec- 
ular systems fl7f| . Third, in this work we have neglected 
decoherence effects (e.g., due to the spontaneous emis- 
sion from the excited state |2)). It would be interest- 
ing to explore how decoherence affects phase control and 
quantum-classical correspondence in this system. 

In conclusion, we have demonstrated optical phase con- 
trol of nonadiabaticity-induced quantum chaos in a A- 
type 3-level system in a two-standing-wave optical lat- 
tice. The results shown in Fig. [| are but samples of 
the observable, controllable, behavior. Further, the func- 



tional dependence on cf> has been exposed analytically, 
which can serve to guide experimental studies of the 
^-dependent regular to chaotic transition. Recent ex- 
perimental progress in atom optics and quantum chaos 
P,p| p^Jl5| , ^8[ | suggests that the results should be experi- 
mentally achievable with existent technology. This work 
was supported by the U.S. Office of Naval Research and 
the Natural Sciences and Engineering Research Council 
of Canada. J.G. is a Henry Croft Postdoctoral Fellow in 
Theoretical Chemical Physics. 
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